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Abstract 

In this letter we compute the exact effective superpotential of N = 1 U{N) supersymmetric 
gauge theories with Nf fundamental flavors and an arbitrary tree-level polynomial superpotential 
for the adjoint Higgs field. We use the matrix model approach in the maximally confining phase. 
When restricted to the case of a tree-level even polynomial superpotential, our computation 
reproduces the known result of the SU(N) theory. 
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In a series of papers [TJ |2] Dijkgraaf and Vafa argued that for a wide class of Af = 1 
U (N) supersymmetric gauge theories the effective superpotential, thought as a function 
of the chiral glueball superfield S, could be computed by means of a simple matrix model 
whose action is the tree-level superpotential of the gauge theory. Their proposal was the 
result of a detailed study of the various dualities between U(N) supersymmetric gauge 
theories, B-model topological strings and matrix models in the large N limit 0J U\ ■ 
This kind of limit was actually an old idea due to 't Hooft jH] who had shown that 
perturbation theory in the large N limit singles out only the planar Feynman diagrams. 
What Diikgraaf and Vafa discovered was that a leading-order perturbative calculation 
via a matrix model could capture completely an exact quantity in the corresponding 
gauge theory, namely the effective superpotential which describes the effects of gaugino 
condensation. In particular they showed that for a theory with a matter field in the 
adjoint representation of the gauge group the effective superpotential is always of the 
form: 

W eff (S)=N^fl (1) 

where J~ x =2 represents the planar free energy from diagrams with the topology of the 
sphere. Then this result was obtained directly using a perturbative field theory approach 
jHj, and through the analysis of the generalized Konishi anomaly [TU| . Many exten- 
sions of this idea have been studied, with the aim to include also quark fields in the 
(anti)fundamental representation [TTJ H21 HB1 EEH UH1 EEB E3 UHj- In particular in [THj it 
was realized that such a generalization could be achieved by taking into account the con- 
tributions from planar diagrams with the topology of the disk (T x= i), i.e. the diagrams 
with one quark-loop boundary [T3|. In this case the relation ([I]) was modified as follows: 

W eff (S) = N d:F x= 2 s {S) + F X=1 (S) (2) 

The validity of this approach was tested considering a theory whose lagrangian contains 
a mass term for the quark fields q (q), a Yukawa interaction with the field in the adjoint 
cf) and a quadratic tree-level superpotential for (p. This model has been studied also for 
different gauge groups [EU EDI EE] and a generalized Yukawa coupling [22] • 

In this letter we focus on a generalization of this kind of models, with a M = 1 U(N) 
gauge theory with Nf flavors and an arbitrary tree-level polynomial superpotential for (ft, 
i.e. 

n+l 

W W = E pk^ k (3) 
fc=l 

We want to obtain the matter contribution JF X=1 (S') for this general case. Following 
we study the maximally-confining phase of the theory and compute T x= \ explicitly as a 
power series in S, i.e. 

F x =i(S) = J2P^ (4) 

j 

We find that the coefficients Pj for given j depend only on the couplings in the potential 
with k < 2j. 
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Our general derivation naturally includes the results obtained by Gomez- Reino [2"31 , 
who solves a SU(N) matrix model using the properties of the factorization of the Seiberg- 
Witten curve. This is a direct consequence of the fact that for a maximally confining 
SU(N) theory only the moduli carrying even indices contribute to the factorization of 
the Seiberg-Witten curve |2Sj- Therefore restricting our result to the case of an even 
polynomial tree-level superpotential we recover the results in [21]. 



We consider a M = 1 U(N) supersymmetric gauge theory obtained by softly breaking 
M = 2 supersymmetry via the introduction of a tree-level superpotential. The action is 
given by 

S S ma n er -\- Sg aU g e -\- S'breafc (5) 



with 



S m atter = j ' d A x J d 2 8 d 2 (e v (f)e Y (f) + qe v q + qe v q 

Sgauge = 27UT J d*X J SB Tr(W a W a ) 
Sbreak = j d * X J d 2 Tr (W tr ee((j), q,<l) 



(6) 



and 



W tree ((p, q,q) = W{4>) + mqq - q<j>q (7) 

where <fi is the matter field in the adjoint representation of U(N), q and q Nf pairs of quark 
fields (with mass m) in the fundamental and anti-fundamental and W a the field-strength 
of the theory. The superpotential W(4>) has the general form given in (j3J). 

At the quantum level the vacua of the theory are determined by the appearance of a 
gaugino condensate described by a chiral superfield 

s = ^ T ^ Q wg (8) 

It is the effective superpotential W e ff(S) which encodes the vacuum structure of the 
theory. 

We follow the matrix model approach of PjJ E] with the generalization to include 
fundamental matter fields 13 , so that we replace <fi with a iV x N hermitian matrix <3> a 6 
(a, 6=1, N), q with a.NxNf matrix Q a a and q with a N f xN matrix Q% (a = 1, Nf) 
and write the matrix integral 



Vol(U(N)) 



J d§ dQ dQ e -^ Wt ^' m (9) 



where Vol{U(N)) is the volume of the gauge group and g s is the string coupling. 

It is well known jTH] that in the 't Hooft limit, where we let iV — > oo, g s —>■ while 
keeping Ng s fixed, the matrix model partition function Z receives contributions from 
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planar diagrams both with the topology of the sphere (x = 2) and the topology of the 
disk (x = 1) corresponding respectively to and 1 quark boundary. Thus, if we call T 
the free energy of the model, we have 

Z = e T « e~^~ 9s (10) 

Now, if we interpret Ng s = S as the glueball chiral superfield of the gauge theory 
(jSj), we are led to the expression in (j2J) The first term in (J2J) is exactly the glueball 
superpotential conjectured by Dijkgraaf-Vafa for theories with fields only in the adjoint 
representation while the second one comes from the extension to include quark fields 
in the (anti)fundamental. 

We want to evaluate this matter contribution T x= \ for an arbitrary polynomial su- 
perpotential for $ of the form (|3*|). As shown in JT] T x= \ can be written in terms of an 
hyperelliptic curve y(x) as follows: 

T x=1 = ~N f [ A \W(x) - y{x)) dx (11) 

where m is the mass of the quarks and Ao a regularization cut-off. In order to integrate the 
hyperelliptic curve y(x) in (fTTj) here we will follow the approach of [231- We study the so- 
called maximally confining phase of the theory so that the hyperelliptic curve degenerates 

as 

y(x) = y/W'(x) 2 - fn-iix) = P n ^(x)^(x-af-^ (12) 

where P n -\(x) is a polynomial of degree n — 1 and a is a parameter that can be shifted 
to in a U (N) theory. The crucial point is to perform the change of variable 

^-f^ + r 1 ) (13) 

in order to expand W(x) and W'(x) in series of £ and 

n+l 

k=l 

n 

W '(x) = W'^it + C'^co + Ys^ + C*) (14) 

fc=i 

In this way one obtains (see [2B] for details) 

n 

y(x) = P n _!(x) v 7 ^ 3 ^ = Ck ^ ~ ^) ( 15 ) 

k=l 

which allows to solve for S 

Ci 

S = I y(x) dx = /i— (16) 
J-ti 2 
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and gives also 



/n+l 
y (x) dx = -ficrfogffl + b + J2b k (£ k 
k=i 



(17) 



Now we need to compute the coefficients bk and Ck in Performing the change of 
variable (pH?j) we obtain 



with the following convention: 



k=l 



j=k 



J 

2m+l 
2 







2m J (2m)\(j - 2m)\ 



(19) 



From ()14|) and ()18|) we learn that: 



n+l 



E 

n+l 

E 



j V27 V V 



(20) 



Following an analogous procedure we obtain the value of C\ which is needed in (II 6|) : 



ci 



2^ 



2 



(21) 



As shown in ED one finds 



T x=x = ~N f l -W(m)+ii Cl log( 



e(A ; 



n+l 



fc=i 



1 + A A 



1 - \ 1 



")' 

to/ 



:22i 



where from ()13|) one can see that: 

Ao 



e(Ao) 

{(to) 



/' 

m 



1 + i/l 



An 



2Ac 
A* 



A r 



to/ 



(23) 
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Using the following identity 



;i + a) '-(i-af = 2yj( 2m * +1 



2m+l 



(24) 



m=0 



we can rewrite 



as 



OA 

- W(m) + 5 logl — ) + b - 2S log 



m 



1 + ^1- (fi/my 



n+l 



k=l r Z=0 



fc 



2i + l 



This equation can be easily expanded in series of /i 2 . Using 

•1 + y/1 - (fi/m) 2 \ ^ 1 1 / 2 / \ / // \ -J 



_ V— — ( 23 V-^-V 

^ 2? 2 2 i I J / \m) 



21 + 1 



1-1^ 

m 



and the expansion in (|20j). finally we obtain 

n+l 



£< 

i>o 



21+1 

yi I 2 



m/ 



2.; 



J 7 y=l 



2 / 



i=i J j'=i 



n+l / n+l 



V— — ( 23 V-^Y 3 ' 

^2j2 2 A j )\m) 



+2 £(£f(f)*(4) £(/ +1 ) B-KTX 



fc=l \ i=fc 



Z=0 



j>0 



(25) 



(26) 



mJ 



(27) 



Our final aim is to write T x =\ as a power series in S. In order to do so we have to work 
on ()27|) showing that it can be drastically simplified. 



Let us consider a matrix model with Nf flavors and a general superpotential for the 
chiral superfield in the adjoint representation of U(N), given by the sum of an even 
polynomial Wi n and an odd one W^n+i : 

n n 

w(x) = w 2n (x) + w 2n+1 (x) = 9 f* 2j + E #T7 x2 ' +1 (28) 

3=1 3 i=i 3 
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We notice that, as explained in the case W(x) = Wi n {x) corresponds to a SU(N) 
theory (for which only the even terms contribute to the glueball superpotential) . The 
SU(N) theory has been considered in [23], where the effective superpotential was com- 
puted order by order in S using factorization properties of the Seiberg-Witten curve. 
Solving the general case in (|2*%|) we will be able to compare our results with those in [23] • 
First we rewrite (|2*7)l separating the even part from the odd one: 



9 -^m* - J2 ^T^ V+1 + S log 



' 2j ^2j + l 

3=1 J 3=1 J 

j=i j \ i=i 

n / n 

+ 2 E( E 



2A r 



+ 



k=l 



J 

2i\2 



U 2j 2* V j 



2j 



m 



+ 



2 ' / 2z 
i — k 



jfe-i 

E 



2=0 



2fc 
2/ + 1 



E< 



22+1 
\ ' / 2 



/i\2(i-fe) 
mJ 



+ 



EE 

fc=0 V i=fc 



g«+i fjVV 2i + 1 

2z + lV2/ V i -A; 



E 

1=0 



2k + 1 
21 + 1 



E<- 

3>0 



21+1 
l I 2 



m 2k-2j+l 2(j-k+i) 



(29) 



In order to obtain .F x =i as a power series of S we proceed in two steps. First we look 
for an expansion in /i 2 . Then we will express /i 2 as a power series in S by means of the 
formulas and (|21j) . To this end we organize eq. (|29jl in the following way: 



^=1 = 



E 



#2j 2 7 



> — -m J — > — 

^ 2 j ^ 2 j + 

3=1 J 3=1 J 



y 9mL m ^i + s log ( 2 M + 

^ 2j + 1 y V m J 



+ e A ^y + E + e c ^ 2 y + E ^(a* 



2 V 



i>2 



i>o 



(30) 



where the seven terms correspond to the seven terms in ([29)1 . 
The A,- coefficients are immediately identified 



Ai = 



92j 1 / 2j 



2j 2V V j 



(31) 



The Bj coefficients are not too difficult to compute and one finds 



3-1 



B i = E 



g 2 i ( 2% \ ( 2(j - i 



i=l 

i<n 



2 2 3 



j -i ) 2(j - i)m 2 0' _i ) 



(32) 
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On the other hand the computation of Cj and Cj is a more difficult task. Let us concen- 



trate on them. In the last two terms in ()29|) we perform the change of index i 
and a shift j — > j — i. In this way we obtain 



i + k 



n j+k 



c ' - EEf 



9%i ( 2i 



1 



\j-i+k 



k=l i=k 

i<n 

n j+k 



2i\i — kJ 2 2i ~ 1 m?^-^ 



k-l 

E 

1=0 



21+1 



2k 

21 + 1 A j -i + k 



s, - EE 



9%i+\ ( 2i + 1 



k=0 i=k 

i<n 



2i + l\ i-k J 2 2i m 2 ^- i )- 1 



\j-i+k 



E 

z=o 



21+1 



2k + 1 
21 + 1 A j - i + k 



(33) 



(34) 



The sum over / in ()33|) and (|3*4*j) can be rewritten in a nice form: 



fc-i 

E 

1=0 



2k 



21+1 



21 + 1 J\j-s + k 



2 2 ( s -i) 



k 



T(k + s-j) 

(j-s + k)\ F(2s-2j + l) 



(35) 



fc-i 

E 

1=0 



21+1 



)-2k-l 



T(2k + 2) T(2j - 2s - 1) 



2k + 1 

21 + 1 A j - s + k ) ~ (j - s + k)\(2k)\ r(k - j + s + l)T(2j - 2s - 2k - 1) 



(36) 



We start computing the lowest terms: 



Co 
C 



ol. k—1 

E ™> 92k \- 
2»-i 2k 2-** 

k=l 1=0 



2k 
21 + 1 



k=l 



E 

k=0 



2 2k 2k 



tE 

1=0 



21 



^ 2k 

n 

^ ;// 2/,- • I 92l.-+\ 



k=0 



2k + 1 



(37) 



In fact these coefficients are exactly cancelled by the first two terms in (|30|) . Moreover we 
want to show that all the other coefficients Cj and Cj, j ^ receive contributions only 
from the first 2j and 2j — 1 terms of the polynomial superpotential respectively. 

Therefore we need to show that every term in (}3*3*j) and (j34j) with an index % > j does 
not contribute. Let us consider i = s for some s > j. It is not hard to see that the 
contribution in (|3~3*)) proportional to gi s can be written as a sum over k = s — j, s 



92s 



E 



2s 



\j-s+k 



k J 2 2s m 2 ^- s ) 



. k=s-j 

In the same way the contribution in 
1 / 2s + 1 \ 1 



fc-i 

E 

1=0 



2k 
21 + 1 



21+1 



j — s + k 



(38) 



proportional to <?2s+i becomes 



92s+l 



E 



2s + 1 



\j-s+k 



k-l 



k )2 2s m 2 ^- s y 



r E 



1=0 



2k + 1 
21 + 1 



21+1 



j — s + k 



(39) 
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Now, setting t = s — j, t > (J5o|) becomes proportional to 

20' + *) \, sk,u , ^ T(A; + 2t) _ 



£i=0 



4^r(j + i/2)r(2t) _ r(2j + i)r(2t + i) _ 
v^./ >■(./) AAA/! 

(2*-* r (i ) r (j + 1/2) - ^r(2j)) s o 



In a similar manner from (}39|) we obtain 

^ T(2A; + 2) T(2j-2s-l) 

22fc Ti - s + 

fc=s-j 



2 2fe (j - s + fe)! r(Jfe - j + s + f ) T(2j - 2s - 2k - l)(2Jfe)! 



which is proportional to 

4^r(j + i)r 2 (j>i)-2r(jofjv / ^(j-2s-2)r(2j)+4-'(i+ s )r(j + i)r(j + i) 



2' 



2' 



We have used repeatedly the following property of the T matrices (for integer j 

1 \J + l l l ) — 22j _! 

This completes our proof. 

At this point we are left with: 

j j-i 

Cj = Cjs92s Cj = Cjs92s+1 

s=l s=l 



where 



c 



AA S+ * / 2s \w 2A; W *±i 



and 



.7* 



fc=l «=o 

A (-)J- 8 + fc / 2s + l \ A / 2Jfe + l W 2I±1 
^ 2 2s (2s + l)m 2 0-«) \ s-k J ^ \ 21 + 1 )\ 



2 2s (2s + l)m 2 0- s ) V s-k J ^A 2/ + 1 / V j — s + A; 

fe=l v 7 1=0 

Using and (J5fijl we find: 



A 



1 /2j 



2 2 J'2jV j 

1 r(2s) / 2{j - s) - 1 

2^-1^ sm 2(i- s) r 2 ( s ) v ./ 1 s " ; 

J r(2a + l) / 2{j - s) - 2 

2 2j-lj s 2 m 2(j-s)-l T 2^ { j- S -l 
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Note that the Aj terms are exactly cancelled by Cjj for j = l,...n computed in ()46p. 
leaving no term linear in /i 2 (then no term linear in S except for the standard piece 
S log(2A /m)). So we can rewrite (}30j) as follows: 



where we have defined 



-N f 
2 } 



i>2 



Cj j >n 

We give explicitly the form of the lowest terms for the coefficients Bf 



B 2 
B 3 
B, = 

and for the coefficients Cj and Cf 
Ci = 

c 2 = 

c 3 = 

c, = 

and 



8 m 2 

3 g 2 3 g 4 
64 m 4 32 m 2 



5 g 2 _^_9a_ _5_J7e_ 
192 m 6 256 m 4 64 m 2 



92 

' 4 



92 



;9i 



16m 2 32 J 

3 g 2 _ 3 g A 

96 m 4 96 m 2 96' 
5 g 2 9 g 4 5 g 6 



:96 



35 



256 m 6 512 m 4 256 m 2 1024' 



C 2 
Ca 







93 

8m 

93 



95 



24m 3 16m 
3 #3 3 # 5 



128 m 5 128 m 3 



5 #7 
128 m 



(47) 



(4f 



(49) 



(50) 



(51) 



As a final step we want to reexpress T x= \ in a power series of S. From (fT6j) and (J2"T|) 
we know that 



i=i 



(52) 
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We want to invert this relation in order to obtain fi 2 in terms of S. We look for an 
expression of the form 



E a mS m (53) 

m>0 

where the coefficients to be determined. Inserting ()53|) into (J52)l we obtain 



j=l m>0 
fn+li 



= ^Ef|(/) E M( j;pi ,p 2 ,...)afar-^ 1+2P2+3P3+ - (54) 

j=l pi,p2,—=0 

pi+P2+---=i 

where we have defined the coefficient: 

j\ 

With a bit of labor one can argue that the general coefficient be expressed recur- 

sively in terms of ai, a,2, a m -\- 
4 

a x = — 

92 



2 / 2(pi + ... + p m _i) \ Pl Pm _! 

a l ••• a m-l 



32 Pl + ••• +Pm-1 

pi+2p2+— +(m-l)pm-i=m 

Pl+P2 + -+Pm-l<[ Ii J i ] 

(pi + ... +p m -i)! 02(pi+...-H» m _i) 



(56) 



(pi)L.(zW-i)! 2 2 (pi+-+p™- 1 ) 
Here we give the expressions of the first coefficients: 

a 2 = -12^7 

#2 

as = +724-40^ 

#2 #2 

a4 = _ 5 4o4 + 600^ - 140^| 

#2 #2 92 

a 5 = +45364 - 7560^ + 2520^ + 120o4 - 504^ (57) 
9 9 2 9% 9l 9 7 2 9% 

Now, using (J55j) and (J5T)Jl we are able to give the expansion of T x= \ in series of S: 



—N f 
2 ' 



S 1 /op ( — - 

V m 



E E(^ +C i +C i) E m 0";pi. 



fc>2 



-,pk 



s k 



pi...p fc =0 

pi+...+p fc =i 
pi+...+fcp fc =fc 
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-Nf 



+ 



—S log ( — - 

2 y V m 



+ 



2g2in 2 



93 

g\m 



S 2 



1 



4 gz 



9i 



+ 77 



2g|m 4 3 (?2 777,3 



5 1 



9% 777,2 



, g3g4 



#5 



V 3 



#3 



+ 18 



6 gfm 6 
9A9b 



4 ^ 



+ 9 



9 9, +12 ^1_ 45 



2 (?2 m2 4 g 4 m 4 

5^6 +20 ^6_ 5 ^ 
2 gjm 



95 



g\m 



4 ^ 

9i m 



g\rn? 



939% 
g\m 



g\m 



g\m 



S 4 + ... 



(5? 



If we extract from (j58j) the even contributions we reproduce exactly what has been com- 
puted in 



In order to have the full expression of the glueball superpotential (j2J) we have to add 
the x = 2 contribution to (|58j) . One can find the implicit solution to this problem in [2E] 
where the x = 2 contribution is given in terms of some parameters which are non-linear 
functions of S. After some algebra we found that the explicit solution can be written as 

N^f = -NSlog{£^) - W S )S> (59) 



j>0 



where 



i~ X (_)m+l m 

W ^J2 -^T E M(m;p!, ...,p j -i)D 2 )1 ...D p j J ~ 1 (60) 



and 



m=l ^ pi,...,pj_i=0 

Pi+...+Pj-i=m 

Pi + ---+0'-l)Pj-i=i-l 



2j [p/2] n p-2q 



<^£f£U)( 7) £ M(p-2 5 ; Pl ,..., B _ 1 ) Z f'...4r 1 ' 



p ^ V 2g / V g . 

p=2 <?=0 pi,...,p.,_i=0 

pi+...+Pj_l=p-2g 



E M(q ] q 1 ,...,q j )D?...Df (61) 



9lv,9j=0 
<?i +...+<?.,■ =g 

(pi+gi)+...+(j-l)(Pi-i+gj-i)+jgj=j 



are defined in terms of some parameters Dj and Zj which can be computed recursively 

k 

D 3 = J2dk Yl MikM,. (62) 

k>0 pi,...,Pj=0 
pi+...+Pj=fc 
pi+...+jpj=j 



11 



and 



7 - ? 93 

Z\ — — Z— g 
92 



2j [p/2] q 
yi 1 ;d=2 F qr=l * * Pi,...,P,=0 



p=2 q=L p 1 ,...,pj~- 

pi + ...+Pj=q 



p-2q+p!+...+jpj 

M(p -2q + Pl + ... + j Pj ; q x , q j _ 1 )z?...zfs 1 1 (63) 

<ji,...,<?j_i=0 
qi+...+q j _ 1 =p-2q+pi+...+jp :j 

ai+---+0'-i)9i-i=i 

The <ij coefficients which appear in (69) and (70) are also given recursively: 

d = -— 
1 293 

2j [p/2] q 

d i s -t E ( 2 P J ( « ) E -.Pi-i) 

y3 P=3 g=0 H H pi,. ..^-1=0 

Pi + ...+Pj-i=<? 
p-2g+pi+...+(i-l)pj-i=i 

(64) 



Explicitly to the fourth order in 5* we obtain 



N 



X=2 



OS 



N 



o 9 ! _ 

3 

95 



3^4 

2<7 2 2 



5- 2 + 



(32 at 

V3 <? 2 6 



-24 



9 04 , , n 9395 

2 92 92 



9l 

10^)53 + (2804 _340^fi + 270Mf!_^4 
3 <? 2 3 / V 3 ^ 2 9 s 2 8 (g 2 ^ 



+ 



+ 



+200^ - 180^^ + 184 

#2 #2 #2 



1 00 M^ + 30^ + 40^1 - 5 4 + .. 

,6 „5 „5 4 q 4; 



</2 



92 



(65) 



Let us notice that the S k term in this expression depends only on the low order cou- 
plings #3, #4, ... , #2fc) in a way parallel to what we have found for the case of fundamental 
matter (see ([58))) . When only an adjoint matter field is present this behavior follows im- 
mediately from the perturbative analysis in ^j. There it was found that contributions to 
the S k term come from a planar diagram with exactly k loops, which can contain vertices 
with at most 2k external legs. 

What we found in (|58|) is the corresponding situation where now the graphs have an extra 
loop of fundamental matter. It is very easy to write down the graphs corresponding to 
every single term in the expressions (|58|) and since the factors of g 2 and m in the 
denominators count the number of adjoint and fundamental propagators respectively. 
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To summarize, we have computed the effective glueball superpotential for a M = 1 
U(N) gauge theory with matter in the fundamental (Nf flavors) and a field in the adjoint 
with a general polynomial tree-level superpotential. This has been done by using the 
technology developed in j2H| for computing the contribution to the matrix model partition 
function coming from diagrams with the disk topology. The full contribution to the 
glueball superpotential is given in the expressions (J58|) and 

Our results generalize the ones obtained in [23] since as discussed in the SU(N) theory 
considered there is equivalent to a U (N) model with an even polynomial superpotential 
for the adjoint field. 

It would be interesting to extend these results to the case of a superpotential admitting 
several vacua. 
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